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Giant enhancement of the thermal Hall conductivity n^y in the superconductor 

YBasCugOy 

Y. Zhangi, N. P. Ong^, P. W. Anderson^ D. A. Bonn^, R. Liang^ and W. N. Hardy^ 
^Joseph Henry Laboratories of Physics, Princeton University, Princeton, New Jersey 08544 
^ Department of Physics, University of British Columbia, Vancouver, Canada. 

(February 1, 2008) 

In high-purity crystals of YBa2Cu307, the quasiparticle (qp) lifetime r and the (weak-field) thermal 
Hall conductivity K^y undergo dramatic increases below 90 K. We present a detailed picture of the 
behavior of Hxy at low temperature, in particular its scaling properties, which are directly relevant 
to the issue of whether Landau quantization of the qp states occurs. 
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The problem of excitations of the superconducting con- 
densate in the cuprates at low temperatures is of strong 
current interest. In a d-wave superconductor, the energy- 
momentum dispersion of quasiparticles near a node is 
Dirac-like. The effect of an intense magnetic field on the 
quasiparticle (qp) states is an interesting open question 
[0-^. Landau quantization of the qp states, first pro- 
posed by Schrieffer and Gorkov jl], has been recently 
re-derived using different arguments (^-^. However, the 
case against Landau-level formation has also been argued 
i|. 

A second problem is the temperature dependence of 
the qp mean- free-path £ (in zero field) close to Tc. Trans- 
port evidence from thermal conductivity Q], microwave 
and teraHertz experiments , and thermal Hall con- 

ductivity |p^-^15|| point to a sharp increase in the qp 
lifetime just below Tc- Recent high- resolution angle- 
resolved photoemission (ARPES) experiments [ [ll|Jl2| 
have started to address the lifetime issue as well, but 
with conflicting results (see below). 

These issues reflect the strong interest in the low- 
lying excitations of the d-wave superconductor. While 
microwave absorption and ARPES experiments provide 
valuable information on the quasiparticles, they are less 
effective in a field. For in-field experiments, teraHertz 
techniques jl^ and the thermal Hall effect |l^-[l5|], in 
particular, have emerged as powerful probes of qp trans- 
port. In a field, the qp heat current develops a transverse 
component that is observed as a thermal Hall conductiv- 
ity Kxy (by contrast, phonons do not display a Hall ef- 
fect since they are charge-neutral). Hence, Kxy selectively 
senses the qp current alone . To fully exploit this tech- 
nique at low temperatures, however, samples with a very 
long £ are needed. 

A recent innovation is the growth, using BaZrCOs 
(BZO) crucibles, of crystals of YBaaCuaOy (YBCO) 
with nearly perfect crystalline order (from x-ray rock- 
ing curves |6|) and very low impurity concentration. 
The step-wise improvement in crystal quality results in 
strong enhancements of the qp lifetime t. A number of 
novel features of qp heat transport become apparent in 
these crystals. The weak-field n^y undergoes a remark- 



able thousand-fold increase between Tc and 30 K. Below 
30 K, the curves of Kxy vs H provide new, specific infor- 
mation on scaling behavior at low T [Q. Both features 
are directly relevant to the two issues mentioned above. 




FIG. 1. (Main Panel) The thermal Hall conductivity 



vs. H in BZO-grown YBa2Cu306.99 (Tc = 89 K) at tempera- 
ture from 85 to 40 K. As T decreases below Tc, the initial slope 
Kxy/B increases sharply. The prominent peak in Kxy below 
55 K is a new feature in BZO-grown YBCO. The inset com- 
pares the zero-field Kxx = Ka in the BZO-grown crystal (solid 
circles) with a detwinned non-BZO grown crystal (open). 

Among the cuprates, 90-K YBa2Cu307 displays the 
largest in-plane thermal conductivity anomaly. In BZO- 
grown crystals, this anomaly is further enhanced, as 
shown in the inset in Fig. The longitudinal ther- 
mal conductivity Kxx (— VT || a) in BZO crystals (solid 
circles) attains a peak value that is ~ 80% larger than 
that seen in typical, non-BZO detwinned crystals (open 
circles). To isolate the qp current, we turn to Kxy The 
main panel of Fi^ |^ displays traces of Kxv vs. field B 
from 85 to 40 K ||l^. As in earlier studies the ini- 
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tial slope n'^y/B = lims^o K^y/B increases very rapidly 
as the temperature T falls below Tc. Further, the curves 
are strongly non-linear in H. Both features reflect a r 
that increases rapidly with decreasing T. An important 
new feature, absent in previous studies, is the prominent 
'overshoot' that produces a maximum in Kxy at the field 
scale Hmax- As T falls below 40 K (see Fig. |), the 
peak continues to narrow. For later reference, we note 
that, over a broad range of temperatures (10 < T <70 
K), Hmax varies as T^. Moreover, at low temperatures 
(T <28 K), the peak magnitude K^y^ also scales as T^. 
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FIG. 2. The thermal Hall conductivity Kxy vs. H in 
BZO-grown YBa2Cu3 06.99 between 35 and 12.5 K. Below 28 
K, the peak value varies as (see text). 

The initial slope n^y/ B, plotted as solid circles in Fig. 
^, undergoes a thousand-fold increase between and 30 
K (the T-linear variation of Kxy above Tc are displayed as 
open circles iQ ) . We now show that this giant enhance- 
ment is driven by a 100-fold increase in the qp lifetime. 

To extract the zero-field mean-free-path (mfp) t from 
K^y/B, we apply the Boltzmann-equation approach p9| , 
which should be valid in the weafc-field regime uj^t <^ 
1 {ujc is the cyclotron frequency). In terms of the 'qp 
heat capacity' Ce = T~^J2k{~(^f /9Ek)E^ where is 
the qp energy, the zevo-H thermal conductivity may be 
written as Ke = Ce{v£)/2, with the group velocity Vk = 
S/Ek/h. (Close to a node k*, the qp energy may be 
approximated as Eq — h^J{vfqiY + (''^A92)^, where Vf 
and wa are velocity parameters normal and parallel to 
the FS, and q = k-k*.) 

The thermal Hall conductivity is related to Kg by n^y — 
Ketanfl. We assume that, in the weak-field limit, the 
thermal Hall angle t&nO is proportional to WcT, viz. 



taii9 = r/ujcT = rii/kpe^, [B ^ 0) 



(1) 



where the magnetic length Ib = y/h/eB. The param- 
eter T] is less than 1 if £ is anisotropic around the FS. 
To obtain tan 6' ||l^, we first fit the profile of Kxx vs. 
H to the empirical expression Kxx{B,T) = k'^(T)/[1 + 
p\B\^] + Ki,g{T), where the background term Ki,g{T) is 
i/-independent and identified with the phonon contri- 
bution. The initial Hall angle is then obtained as |l^] 
tan6' = liiJiB^o Hxy{B)/[KxxiB) — Kbg]- This procedure 
allows us to extract tan^ (hence, i using Eq. |l|). 

As a consistency check, we adopt a second way to ob- 
tain £ from kI,, that relies on measurements of the elec- 
tronic heat capacity Cg. Using Eq. 0^, we may write 



-xy 



(2) 



10" 



^ 10- 



> 10- 



10" 

















{ 












,1 , , , ' 



10 



100 



r(K) 



FIG. 3. The T dependence of the initial Hall slope K^y/B 
in BZO-grown YBCO (solid circles). Between and 30 K, 
K-xy/B increases by 10^. The l/T dependence of K^y/B above 
Tc (measured in a non-BZO grown YCBO) is shown as open 
circles. The inset shows a qp energy contour on the Dirac 
cone. Group velocities on the particle- (p) and hole-like (ft) 
branches are indicated. 

In a d-wave superconductor, Ce = olcT^ for T < Tc. 
Using the measured value ac — 0.064 mJK^'^mol^^ [po[ , 
we may invert Eq. ^ to find £. We find that the values 
of £ obtained from the two methods share the same T 
dependence, but differ by a fixed factor of 1.5 if 77 = 1. 
By adjusting t] to 0.6, we obtain numerical agreement 
between the two methods. 

Figure ^ shows the T dependence of £ derived from the 
two methods. The agreement between the two sets of 
data is evidence that our assumption Eq. |l| is physically 
reasonable. Remarkably, between Tc and 20 K, the mfp 
increases by a factor of ^120 from 80 A to 1 micron. In 
the expanded scale, we show that this increase is abrupt. 



2 



starting slightly below 
1.0 



a 




01 



r(K) 

FIG. 4. The zero-field mean-free-path i extracted from the 
weak-field Hall angle tan 6' (open circles), and from Eq. |^ 
(closed). The equivalent values of 9/B are shown on the 
right scale. The symbols (x) represent tan6' measured in 
a non-BZO detwinned YBCO crystal (Krishana et al. [14]). 
The expanded scale (dashed lines) highlights the steep in- 
crease below Tc- To extract I, we used the values r] = 0.60, 
Vf = 1.78 X 10^ cm/s, and kf = 0.8A'^. 

[For comparison, taiiO measured previously in a non- 
BZO crystal [|l5[ is shown as x . Based on the higher 
sensitivity and broader range in T in the present exper- 
iment, we now conclude that tan^ does not lie on the 
extrapolated curve for the electrical Hall angle tan^e-] 

Beyond the weak-field regime, we need a fully micro- 
scopic description of the qp thermal Hall current to prop- 
erly analyze n^y vs. H. As the theoretical situation 
is unsettled, we adopt instead scahng arguments [^ . 
This approach reveals some rather striking features in 
the data. 

For states close to the node k*, the linear energy dis- 
persion E — hvq {v is an average velocity) implies a 
general relation between ksT and the magnetic length 
Ib Sit a, characteristic field scale Bg, viz. 



kgT — hv 



h 



(3) 



In addition to this general relation, Simon and Lee 
^ have proposed that, at low T (<30 K for YBCO), 
the magnitude of K^y should scale as 



K,yiH,T)^T^F,yiVH/aT), 



(4) 



where a = ks/vVeh, and Fxy{u) is a scaling function of 
the dimensionless parameter u — ^/H /aT . Hence, plots 



of Kxy/T"^ versus \Jli jT should collapse to the universal 
curve Fxy{u). 




+ 


15K 


■ 


17.5 





20 


A 


22.5 


V 


25 


— • — 


27.5 


— A — 


30 


— T — 


35 



0.12 



(^„/^)"'/r(T"'/K) 



FIG. 
(Eq. [ 
versal' 



Simon-Lee scaling plot of K^y/T^ versus y/H/T 
^. Below 28 K, the curves collapse onto a 'uni- 
curve F^yiu). Above 28 K, scaling is violated. 
However, the peaks still occur at the same a;-coordinate 
{VH^/T = 0.042). 

Hare- 



The arrows indicate the field scale 



We proceed to plot our results in this way in Fig. 
While the curves above 28 K are spread out, the ones 
below collapse onto a common curve for H < Hmax ■ The 
data taken at 25 K (and below) collectively determine the 
form oi Fxy{u). Its most notable feature is the nominally 
straight segment that extends from it ~ to just below 
Wo = ^JUmaxloT, i-e. F^yiu) u for < u < uq. This 
simple form for F^y implies that, below 25 K and for 
Fl < Hmax-, Hxy rcduccs to the form 



,iH,T)^CoTVH, 



(5) 



where the constant Cn — 1.51 x 10 ^ in SI units. Re- 
markably, when Eq. |5| applies, the magnitude of Kxy is 
just proportional to tVH, and is insensitive to all trans- 
port quantities such as £ and 6. This interesting result 
has not been anticipated theoretically. 

At larger values of u, F^y attains a maximum value F^y 
before falling slowly. The dependence of the peak 
value ti^y^ noted earlier in Fig. |[ is now seen to be 
a simple consequence of scaling behavior (i.e. k" 

Above 28 K, Simon-Lee scaling no longer holds. Three 
field regimes are now apparent. In weak fields {0 < H < 
Hx), Hxy is strictly linear in H. Above H^, we enter a 
regime reminiscent of the y/H behavior at low-T (the H- 
linear regime is too small to resolve below 28 K). This 



^max 
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intermediate regime appears as straight-line segments in 



xy occurs at the same x-coordinate in Fig. ||, i.e. 

ynax for Bs in Eq. |, 



Fig. |5|. Finally, closer to Hmax, Kxy deviates from V-ff 
behavior, and goes through a broad maximum. Surpris- 
ingly, as noted earlier, the weaker scaling relation in Eq. 
^ continues to hold: Between 15 and 70 K, the maximum 
in K 

y/H^ax = 0.042 T. Substituting H, 
we find that v ~ 8.0 x 10^ cm/s, which is close to the 
geometric-mean velocity y/VfVA 6.8 x 10^ cm/s (with 
Vf = 1.78 X lO^cm/s Q and Vf/vA ~ 7). 

We may estimate semiclassically the time that an ex- 
citation, in a field, takes to move from 1 to 2 along the 
arc (inset. Fig. ||) by = {h/eH)J^ ds^ \ Vk |~^, where 
Sk is the arc-length. For this time to equal r, the field 
required is Hare = ttE /{evAVfr). Using the measured 
£ ~ vfT at each T and setting E = ksT, we indicate 
Hare as arrows in Fig. |[ This rough estimate shows 
that the peak is related to the maximum arc length of 
the dominant energy contour on the Dirac cone. Hence, 
a detailed analysis of the Hall results should shed im- 
portant light on the current debate about how vortices 
affect the qp spectrum. The presence of Landau Levels 
or absence will presumably have a large effect 
on Kxy. Moreover, the direct measurement of the scaling 
function F^y (Fig. |^) together with the other scaling fea- 
tures uncovered should stringently narrow the range of 
possibilities in this interesting problem. 

The new results on K^y also bear on the issue of the 
change in qp lifetime at Tc. As discussed, £ derived 
from transport undergoes a steep increase just below 
|f7|-p|,p^. Recently, ARPES has attained enough resolu- 
tion to probe the qp spectral peak along the nodal di- 
rection in Bi2Sr2CaCu208. Valla et al. |l^ find that the 
width Afc I/Iarpes) retains its T-linear dependence 
across Tc (near Tc £arpes — 25 — 30 A). This appears to 
be in striking contrast with the transport results. How- 
ever, Kaminski et al. resolve a new feature of the qp 
peak that suggests a rather different picture. Below Tc, 
the qp width becomes dramatically narrower provided it 
is probed within 60 mV of the Fermi energy. They infer 
that well-defined qp states at the nodes exist only below 
Tc. The steep increase in £ shown in Fig. |4| is in close 
agreement with Kaminsky et al. The data in Fig. ^ show 
that £ increases to ~ 1/im below 20 K (implying a peak 
200 times narrower than the peaks resolved in the current 
ARPES studies). Hence, in high-purity YBCO, there are 
exceedingly sharp qp peaks in the spectral function that 
remain to be resolved and investigated. Understanding 
the abrupt appearance of the qp state below Tc, as im- 
plied by the steep increase in £ and K'^y/B near Tc, seems 
a key problem in the cuprates. 
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